Conductance of the Single Electron Transistor for Arbitrary Tunneling Strength 
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We study the temperature and gate voltage dependence of the conductance of the single electron 
transistor focusing on highly conducting devices. Electron tunneling is treated nonperturbatively 
by means of path integral Monte Carlo techniques and the conductance is determined from the 
Kubo formula. A regularized singular value decomposition scheme is employed to calculate the 
conductance from imaginary time simulation data. Our findings are shown to bridge between avail- 
able analytical results in the semiclassical and perturbative limits and are found to explain recent 
experimental results in a regime not accessible by earlier methods. 
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The single electron transistor (SET) has become a pri- 
mary device for applications of the Coulomb blockade 
effect . The modulation of the current in dependence 
on the gate voltage has been used to build highly sen- 
sitive electrometers and thermometers and can serve as 
switch or amplifier. Since the relative current modulation 
is most pronounced at small transport voltages where 
also heating effects become negligible, the quantity of 
chief interest is the linear conductance. There are two 
relevant dimensionless parameters Q; the dimensionless 
parallel conductance g of the two tunnel junctions mea- 
sured in units of the conductance quantum Gk = e^//i, 
and the dimensionless inverse temperature l3Ec relating 
(3 = l/ksT with the charging energy Ec- While re- 
cent experiments have explored a large regime of 
tunneling conductances and temperatures, available the- 
oretical work is restricted to limiting cases where one of 
the parameters g or PEc is small. Approaches for the 
weakly conducting regime g~\ [^,|| start in the charge 
basis and treat tunneling as a weak perturbation, while 
the semiclassical approach starts from the canon- 
ically conjugate phase representation and treats small 
fluctuations about the classical paths which is adequate 
for f3Ec~l- So far no theoretical predictions are avail- 
able for the regime of highly conducting systems at low 
temperatures which is within reach of state-of-the-art 
experiments. 

In this work we study the linear conductance of the 
SET for arbitrary parameters. Electron tunneling is 
treated nonperturbatively by means of a path inte- 
gral Quantum Monte Carlo (QMC) approach employed 
to evaluate the current-current correlation function in 
imaginary time. We present an improved singular value 
decomposition (SVD) scheme utilized to determine the 
conductance via an inverse Laplace transformation. This 
allows us to obtain results for the entire range of param- 
eters of experimental interest. 

Specifically we consider a SET consisting of two tunnel 
junctions with conductance Gi, G2 and capacitance Ci, 



C2, respectively. The junctions in series are biased by a 
voltage source V and the island between the tunneling 
barriers is coupled to a gate voltage Ug via a capacitance 
Cg. From the Kubo formula, we write the dc conduc- 
tance in terms of the current-current correlation function 
F{t) = (/(t)/(0)) in imaginary time r as M 
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Once one has calculated the correlator F{t) analytically, 
it is a rather easy task to continue its Fourier transform 
to get the frequency dependent response. However, when 
numerical techniques are used to evaluate F{t), the con- 
tinuation becomes a serious problem: The Fourier trans- 
form is known only for a discrete set of equidistant Mat- 
subara frequencies Vn = 2Trn/h(3. Inserting the spectral 
representation of the correlation function H 



dLoF{uj)e 



with spectral density F{uj) > into 
conductance 
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we find for the 



G=|F(0). 



(3) 



To relate this with the numerically determined current- 
current correlation function F(t), one has to invert the 
two-sided Laplace transformation (^. This operation, 
however, is known as an ill-posed problem: small errors 
in F{t) may cause large deviations in F{uj). One has to 
use a regularization scheme to get rid of highly oscilla- 
tory functions. Here, we employ SVD ||l^ regularized by 
constraints explained below. Before addressing this, we 
present the QMC algorithm. 

Explicitly, the current correlator may be written [0 



F{t) = 4^Gcia(r)(cos[(^(T) - (^(0)]) 



(4) 
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where Gci = G\G2l{G\ + G2) is the classical high tem- 
perature conductance. a(r) describes electron-hole pair 
propagation for electrons and holes created in different 
electrodes, and may be represented in the form (0) 



air) 
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with the spectral fimction 



27r 1 



Finally, the average 

A(t) = (cOs[(y9(T) 



(5) 
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describes coupling between electron-hole pair excitations 
and the electromagnetic degrees of freedom, where the 
time derivative of the phase field ^pir) is related to the 
voltage across the SET. To calculate A(t), we use a for- 
mally exact path integral representation Q 

A{t) = i 1 i?[^]e-*^M cos[vp(t) - (^(0)] , (8) 

where Z is the partition function and SYp\ the Euclidian 
action that splits into S\if\ — ScYp\ + •S'tI'/']- The first 
term 



dr 
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where Ec — e^/2C is the charging energy with the is- 
land capacitance C = Ci -\- C2 + Cg, describes Coulomb 
charging in presence of a gate voltage Ug = eUgjCg. The 
second contribution 



ri/3 



St{v\ = 2.g / drdr'air - t') sin^ 



ip{t) - ipir') 
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is due to electron tunneling where g = (Gi -I- G2)/Gk is 
the dimensionless parallel conductance of the SET. 

To implement the topology of the configuration space, 
we introduce winding numbers k labeling different classes 
of paths running from (^(0) = to (p{fif3) = 27rfc. Chang- 
ing the variables, (^(t) = i^kT + (^{T) with ^(0) = ({hf3) — 
0, one may integrate out the gate voltage dependent term 
in Eq. (^, and write 
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with the unnormalized average over a positive measure 



(W) = 
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Here, S'o[C] is given by Eqs. (g) and (|10|) with the gate 
voltage set to zero. Since 5'o[C] is real, we may employ 



standard Metropolis sampling techniques to evaluate the 
path integral numerically. For the temperatures explored 
the data were found to be saturated at Trotter number 
N « 5pEc- For higher temperatures we have chosen 
N > 100 to keep the resolution of the discretized cor- 
relation function high enough for the Laplace inversion. 
The QMC data show statistical errors of at most 3% for 
all temperatures and gate voltages. The denominator 
((exp(— 27rzfcng))) in Eq. ( pT| ) is an averaged sign factor 
that always exceeds 0.04 for the results shown below. 

To analyse the data we start by noting that A{t) also 
has a representation of the form (|^). In terms of the 
spectral functions we then find 



Fiu) = 2Gc 



duj'a{uj - uj')A{uj') . 



(13) 



In view of the symmetry A{t) = A{h(3 — r) equivalent to 



detailed balance ^(—1^) 



A(uj) , we may introduce 



the positive and symmetric spectral function 



A'{uj) ^ A{uj) 
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and write 



A{t) = — / dwA'{uj)Kp{oj,T) 



with the integration kernel 



(14) 
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Now, we may restrict the time domain to r G [0,h(3/2]. 
To invert the integral transform ( p^ ) consider the opera- 
tor )Cp with kernel Kf}{uj, r) 



gir) = (/C^g)(r) = / dcjKp{LU,T)g{u;) , (17) 
Jo 



and it's adjoint K-g 
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Clearly, the operation ICplCp is self-adjoint and positive. 
The eigenvalue equation 



(19) 



can be solved numerically by discretizing the correspond- 
ing integral equation. The resulting eigenfunctions form 
an orthonormal basis {/n(T)} in L^[0, ?i/3/2]. Now, the 
spectral function A' [uj) can be expanded as 



(20) 



where the 
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are expansion coefBcients. Eq. ( poD explicitly represents 
A{u!) in terms of A{t) and serves as inverse of the integral 
transform ([T5|). Since with increasing n the eigenvalues 
/i„ vanish overexponentially, for large n highly oscillatory 
parts become more pronounced when the expansion co- 
efficients c„ do not vanish fast enough. This exhibits the 
origin of the ill-posed problem: a small error becomes 
amplified by a large inverse singular value l//i^ and thus 
a regularization of the procedure is required. 

First, the integral in Eq. ( p^ ) is cut off at a frequency 
Wmax which is chosen so large that the singular functions 
do not alter appreciably with this cutoff. Second, the 
sum in Eq. (pOl) is truncated at rimax- To fix the value of 



'^max we use the noise level 
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of the data where N is the Trotter number and cr(Ti) 
the standard deviation of the data A{Ti). This has to be 
compared with the error stemming from the truncated 
expansion in orthonormal functions 
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To ensure that we do not include fine structure origi- 
nating from the noise of the QMC data, we fix rimax by 
demanding e„„^^ « cr. 

An essential step is to use additional information to 
regularize the SVD scheme. We know that A' is positive, 
symmetric, vanishes for large frequencies, and must obey 
the sum-rule 



dujA'{uj)- 
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= A(0) = 1, 
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To incorporate this, we make the ansatz 

M 
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where M > nmax- The coefficients a„ are determined 
from minimizing a functional according to three con- 
straints. The proximity of the solution to the QMC- 
correlation function is obtained by minimizing 



prox 



fii, 



,an„„x) = E 
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where the c„ are the expansion coefficients (^^ and 
where the points are weighted according to errors in the 



QMC data. To achieve positivity we minimize the inte- 
gral 



$pos(ai, ...,aM) = / dLjA'j^j{uj)'^w{Lj) , 



(27) 



where the weight function vanishes in the region where 
the result ( po[ ) is positive. To account for the sum-rule 
we minimize the functional 



$sum(ai, um) 



1 



du^^^'^'tW . (28) 
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Overall we have to consider the linear combination 



$(ai, ...,aM) = * 
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with two regularization parameters A and fj,. Differen- 
tiating ( p9|) with respect to a„ we get a system of lin- 
ear equations that may be solved. We have tested the 
method with analytical functions corrupted by Gaussian 
noise to optimize the regularization parameters. 

After extracting A'{lu) from the QMC data via SVD, 
wc obtain F{u}) from the convolution ( p^ and the con- 
ductance from Eq. (||). Varying the gate voltage Ug 
the conductance shows a periodic behavior and we may 
define a maximum Gmax = GL _i and a minimum 

I 'tg — 2 

— G\ng=o conductance in dependence on tempera- 
ture. 

To estimate error bars we add a further functional to 
(p9|): We try to minimize the distance of F(ijj) to a fictive 
either too low or too high value F^{u!). This fictive value 
of the conductance still has to be in accordance with the 
other three constraints and leads to an upper and lower 
bound of possible conductance pT| |. The error bars for 
the conductance found by this procedure lie all below 3%. 
Since they are of the order of the symbol sizes they will 
be omitted in the figures. 

As a first test of the approach we compare our nu- 
merical data with second order perturbative results 
and with semiclassical findings j? ,P . The comparison in 
Fig. |l| shows that our results agree well with both lim- 
iting analytical theories in their ranges of applicability 
and therefore bridge between the regions of parameters 
covered by analytical work. 

Recently, P. Joyez et al. j^] and D. Chouvaev et 
al. Q have measured the conductance of SETs in the 
moderately-to-strong tunneling regime. The strong tun- 
neling data by Chouvaev et al. cover only a restricted 
range of temperatures and were found to be in agreement 
with semiclassical theory. In contrast, some of the data 
obtained by Joyez et al. are not only outside of the range 
of the perturbative theory but also explore temperatures 
low enough to show strong deviations form semiclassical 
results. 

However, the experimental data in Ref. Q for strong 
tunneling "ij = 7.3" are found to deviate significantly 
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from our predictions for g = 7.3 (data not shown). There- 
fore, we have reexamined the parameters given in Ref. : 
The bare charging energy Ec was determined indepen- 
dently from resonances in the superconducting state that 
should be reliable within errors of about 10%. The series 
conductance was measured at sufficiently high tempera- 
tures {l3Ec — 0.18) fixing Gci up to a small error. Assum- 
ing identical tunnel junctions, the dimensionless parallel 
conductance was estimated as g — AGc\/Gk- This as- 
sumption was consistent with a fit (dashed line in inset 
of Fig. 1^) of the high temperature data to the prediction 
G/Gci = l-f3E*(j/'i with 



E*c/Ec - 1 - (9gC(3)/2^^ + c)(3Ec 



(30) 



where c — 0. Yet, the correct result reads c — 1/5 pjl2| 
leading to the dashed-dotted line in the inset which is 
inconsistent with the data. Even for the highest temper- 
atures explored experimentally, the data should be com- 
pared with the full semiclassical prediction from Ref. . 
As can be seen from the inset in Fig. ^ the high tempera- 
ture data are only consistent with a parallel conductance 
near g = 10 implying an asymmetry of the transistor. 
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FIG. 1. Predictions for the maximum and minimum linear 
conductance of a SET for dimensionless parallel conductance 
g = 2.5 compared to second order perturbative results Q and 
for g — 25 compared to semiclassical findings 

In Fig. H we show for g — 10 QMC predictions for max- 
imum and minimum linear conductance in dependence 
on the inverse temperature. They are compared with the 
experimental findings from Ref. for "g = 7.3" and are 
found to be in good agreement. The error bounds of the 
numerical data are again within the symbol size. The pa- 
rameters for the QMC were not adjusted to improve the 
global fit but coincide with those extracted from the high 
temperature analysis. Remaining deviations between the 
experimental data and our predictions may result from 
the uncertainty of the parameters. 

In summary we have presented a general approach to 
calculate the conductance of the SET for arbitrary pa- 
rameters. We found good agreement with limiting ana- 
lytical theories and could explain experimental findings 
for highly transmitting tunnel junctions. As one of the 



results of this study future experiments in the strong tun- 
neling regime should explore the high temperature be- 
havior of the SET more carefully to allow for a unique 
parameter identification within the range of validity of 
the semiclassical theory. 
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l//3Ec PEc 
FIG. 2. Predictions for the maximum and minimum linear 
conductance of a SET for dimensionless parallel conductance 
g — 2.5 and g = 10 compared with experimental findings from 
Ref. Q. In the inset the experimental data of the renormal- 
ized charging energy Eq/Ec for = 7.3" are shown and 
compared with high temperature approximations (see main 
text) . 
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